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Intraband and interband contributions to the current-driven spin-orbit torque in magnetic ma¬ 
terials lacking inversion symmetry are theoretically studied using Kubo formula. In addition 
to the current-driven field-like torque Tfl = tflhi x u so (u so being a unit vector determined 
by the symmetry of the spin-orbit coupling), we explore the intrinsic contribution arising from 
impurity-independent interband transitions and producing an anti-damping-like torque of the form 
Tdl = tdlhi x (u so x m). Analytical expressions are obtained in the model case of a magnetic 
Rashba two-dimensional electron gas, while numerical calculations have been performed on a di¬ 
lute magnetic semiconductor (Ga,Mn)As modeled by the Kohn-Luttinger Hamiltonian exchanged 
coupled to the Mn moments. Parametric dependences of the different torque components and sim¬ 
ilarities to the analytical results of the Rashba two-dimensional electron gas in the weak disorder 
limit are described. 

PACS numbers: 72.25.Dc,72.20.My,75.50.Pp 

I. INTRODUCTION 

Magnetization dynamics driven electrically by spin- 
polarized currents through the spin transfer torqueir— has 
attracted considerable attention due to its applications in 
memory and logic spintronic devices —£ An alternative 
mechanism, the spin-orbit torque (SOT), has been re¬ 
cently proposed as a means to control the magnetization 
of a single ferromagnetic^— or even antiferromagnetio^ 
layer without the need of an external spin-polarizer. The 
SOT arises from the interaction between the nonequi¬ 
librium spin density of carriers and the local magneti¬ 
zation. The non-equilibrium spin density results from 
the transfer of angular momentum between the spin and 
orbital degrees of freedom of the carriers— r— The 
SOT requires magnetic structures with strong spin-orbit 
coupling and inversion symmetry breaking. Initially ob¬ 
served in epilayers of (Ga,Mn)As dilute magnetic semi¬ 
conductors (DMSs) with bulk inversion asymmetry in 
their strained zinc-blende crystal^ - — this effect was 
soon widely confirmed in metallic bilayers with struc¬ 
tural inversion symmetry breaking—In general, the 
SOT observed experimentally possesses two components, 
a field-like torque Tfl = tflih x u so odd in the magne¬ 
tization direction m and an anti-damping-like torque^! 

Tdl = TDLm x (u so x m) even in m. Here, u so is a 
unit vector determined by the symmetry of the struc¬ 
ture and the current direction,— and tfl and tdl are 
the magnitudes of the field-like and anti-damping-like 
torque, respectively. These torques are also commonly re¬ 


ferred to as the out-of-plane and in-plane torques, respec¬ 
tively, with respect to the (m, u so ) plane. The direction 
of the field-like (out-of-plane) and anti-damping-like (in¬ 
plane) torques and their detailed angular dependence^ 
depend on the crystal structure, while their magnitude 
has been shown to strongly depend on the materials 
considered —~— 

Two main mechanisms have been invoked to ex¬ 
plain the origin of the current-driven torques in non- 
centrosymmetric ferromagnets. In the first scenario, 
the lack of inversion symmetry enables the inverse 
spin galvanic effect— 2 (ISGE), i.e. flowing current di¬ 
rectly produces a nonequilibrium spin density <5S lo¬ 
cally, whose direction is determined by the symmetry 
of the spin-orbit coupling. Recently, it has been pro¬ 
posed that in non-centrosymmetric magnetic materials 
this nonequilibrium spin density may exert a torque on 
the magnetization^— iiS T = ( 2,/ ex / hjN rn ) m x <5S. Here, 
7 is the gyromagnetic ratio, N m the density of mag¬ 
netic moments and J ex the exchange coupling (having 
the dimension of energy) between the itinerant electron 
spins and the local magnetization M = M s m which, 
in this article, is assumed to arise solely from localized 
magnetic moments p so that the saturated magnetiza¬ 
tion M s = pN rn . This is the essence of the ISGE- 
induced SOT. Alternatively, in ferromagnets adjacent 
to a heavy metal, it has also been proposed that the 
spin Hall effect (SHE) present in the heavy metal may 
inject a spin-polarized current into the adjacent ferro- 
magnet, exerting a spin-transfer torque (STT) on the 
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magnetization— i 21 i 30 

A current debate aims at identifying the interplay be¬ 
tween these different mechanisms and their impact in 
terms of current-driven spin torque. In the simplest 
physical picture, SHE induces an anti-damping-like STT, 
while the SOT reduces to a field-like torque generated by 
ISGE— However, it has been recently proposed that the 
incomplete absorption of the SHE-induced spin current 
by the ferromagnet (or, equivalently, the non-vanishing 
imaginary part of the interfacial spin mixing conduc¬ 
tance) may result in a field-like STT component— Sim¬ 
ilarly, in the context of ISGE-induced SOT, recent the¬ 
ories have suggested that spin relaxation and dephasing 
may also lead to a correction in the SOT in the form 
of a anti-damping-like component— -1 - j n R e fs. (T2| and 
ITgil . the anti-damping-like SOT term arises from the elec¬ 
tron scattering-induced spin relaxation. In Ref. [30l the 
semiclassical diffusion formalism was used, whereas in 
Refs. ITT] and [3, the anti-damping-like SOT is obtained 
within a quantum kinetic formalism. It is ascribed to 
spin-dependent carrier lifetimes^ or to a term arising 
from the weak-diffusion limit, which in the leading order 
is proportional to a constant carrier lifetime^ 14 

Intriguing material-dependence of the SOTs has been 
unravelled in various experiments keeping the debate on 
the origin of these components open. 2 - - — The difficulty 
in determining the physical origin of the torques partly 
lies in the complexity of the ultrathin bilayer consid¬ 
ered, involving both bulk and interfacial transport in the 
current-in-plane configuration. First principle calcula¬ 
tions have indeed pointed out the significant sensitivity 
of the torques to the nature of the interfaces— 

In a recent publication, Kurebayashi et al— investi¬ 
gated the SOT in a bulk DMS. They observed a large 
anti-damping-like torque that is not ascribed to the 
SHE since no adjacent spin-orbit coupled paramagnet is 
present. It was then proposed that such a torque has 
a scattering-independent origin in the Berry curvature 
of the band structure, in a similar spirit as the intrinsic 
SHE was introduced about ten years agOi 33 i 34 

In this paper, we present a systematic theoretical study 
of SOTs arising from the ISGE and Berry curvature 
mechanisms in a spin-independent relaxation time ap¬ 
proximation. We focus our attention on the current- 
driven spin-orbit field (called the SOT field), h so , pro¬ 
ducing the spin-orbit torque T = M x h so . These 
SOT fields have an in-plane component of the ISGE 
origin^ hjj° = t F lu so [i.e. lying in the (m, u so ) plane 
and producing an out-of-plane torque] and also an in¬ 
trinsic contribution arising from interband transitions. 
The latter— produces an out-of-plane field of the form 
h™ = TDLU S oXm [i.e. lying perpendicular to the (m, u so ) 
plane]. Analytical expressions are obtained in the model 
case of a magnetic Rashba two-dimensional electron gas 
(2DEG), while numerical calculations are performed on 
DMSs described by the kinetic-exchange Kohn-Luttinger 
Hamiltonian— Parametric dependences of the different 
torque components and similarities to the analytical re¬ 


sults of the Rashba two-dimensional electron gas in the 
weak disorder limit are described. 


II. NON-EQUILIBRIUM SPIN DENSITY: 

INTRABAND AND INTERBAND 
CONTRIBUTIONS IN KUBO FORMULA 

In the present study, we start from a general single¬ 
particle Hamiltonian 

H sys = Hq + Hsoc + Hex + Vjmp(r) — eE • f, (1) 

where the first term includes the spin-independent ki¬ 
netic and potential energies of the particle, the second 
term denotes the coupling between the carrier spin and 
its orbital angular momentum and the third one repre¬ 
sents the interaction between the spin of the carrier and 
the magnetization of the ferromagnetic system. Below, 
we refer to these first three terms as to the unperturbed 
part of the Hamiltonian. The fourth term is the impurity 
potential and the fifth term is the electric field applied 
through the system. Impurities are treated within the 
constant relaxation time approximation while the elec¬ 
tric field is treated within the framework of the linear 
response theory. As discussed below, this electric field 
has two distinct effects on the electronic system: (i) it 
modifies the carrier distribution function from its equilib¬ 
rium Fermi-Dirac form and (ii) it distorts the carrier wave 
functions. The former leads to intraband ISGE contri¬ 
butions, while the latter is responsible for the interband 
(Berry curvature) contribution. To calculate the SOT 
field, we evaluate first the nonequilibrium spin density 
<5S using the Kubo formula 

SS = 2^ Re ^ ^ ka l S l^ k ^^’ kfe l E ' V’ka) 

k ,a,b 

(2) 

where G^ a = (G^J* = 1 /{E F - E ka + *T), E F is the 
Fermi energy, E^ a is the energy dispersion of band a, V 
is the system volume and T is the spectral broadening due 
to the finite lifetime of the particle in the presence of im¬ 
purities. The Bloch state I'f/’ko) in band a can be found by 
diagonalizing the unperturbed part of the Hamiltonian in 
Eq. ©. This expression contains both intraband (a = b ) 
and interband (a b ) contributions to the nonequilib¬ 
rium spin density. Numerical results in Section [IV Bl are 
calculated with the above equation. 

In order to understand the numerical results, Eq. © 
can be rewritten^ as <5S = <5S intra + <5S 1 1 nter + <5Si, nter 
when weak impurity scattering (namely, small spectral 
broadening, T —>■ 0) is assumed. The three contributions 
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are 


SS lntra = V’ka) 

k.a 


(5Si nter 


<5S 1 2 nter 


x S(E ka - E f ), (3) 

~yYl 2Re KV’afcl®IV’6fc){V , 6fc|E ■ vIV’afe)] 

k ,a^b 


T(E ka - E^ b ) 

[(Eka -E kh ) 2 +r 2 ] 2 /ka 


/kb)- 


eh 

V 


^ Im[(^ka|s|V , kb)(^kb|E • vIV'ka)] 

k ,a^b 


( 4 ) 


X 


r 2 - (E ka - Em ) 2 

l(E ka - E kb ) 2 + r 2 ] 2[Jka Jkb> 


( 5 ) 


The first term, Eq. ©, is the intraband (a = b ) contri¬ 
bution arising from the perturbation of the carrier distri¬ 
bution function by the electric field. It is proportional to 
the momentum scattering time (r = h/2T) and is there¬ 
fore an extrinsic contribution to the nonequilibrium spin 
density (i.e. it is impurity-dependent). The second and 
third terms are interband (a ^ b ) contributions arising 
from the perturbation of the carrier wave functions by 
the electric field. The second term, Eq. a, is inversely 
proportional to the scattering time, i.e. it vanishes in the 
clean limit. The third term, Eq. is independent of the 
scattering in the weak scattering limit {E ka — E k b T), 
i.e. is an intrinsic contribution to the nonequilibrium 
spin density. The formalism described above is the estab¬ 
lished linear response theory of a translationally invari¬ 
ant system and has been exploited, for instance, in the 
context of the spin Hall^S and anomalous Hall effect— 
Nevertheless, the distinction between these different con¬ 
tributions is particularly important in the case of the 
SOT since these terms give rise to different symmetries 
of the torque. 

The concept of intrinsic SOT is illustrated in Fig. [T] (see 
also the discussion in Ref. [32). FigureQJa) represents the 
Fermi surface of a non-magnetic Rashba 2DEG under the 
application of an external electric field E. At equilibrium 
(E = 0) the spin direction (pink arrows) is tangential to 
the Fermi surface (grey circle) at all k-points, and the 
total spin density vanishes. Applying the electric field 
accelerates the electrons on the Fermi surface and they 
feel a modified spin-orbit field <5B oc z x p = — z x eE 
(thick blue arrow) around which the spin momenta (red 
arrows) start to precess. For a non-magnetic 2DEG, the 
resulting spin density vanishes while a non-zero trans¬ 
verse spin-current is generated by this mechanism. This 
is the origin of the intrinsic SHE. 

Now, let us consider the case of a magnetic Rashba 
2DEG in the strong ferromagnetic limit [Fig. [ljb)]. At 
equilibrium, the spin momenta (pink arrows) are approxi¬ 
mately aligned along the magnetization direction m (thin 
black arrow) for all k-points of the Fermi surface (grey 
circle). Under the application of an external electric 
field, the spin momenta (red arrows) precess around <fB 


(thick blue arrow) resulting in a non-vanishing spin den¬ 
sity. Following the convention adopted in Fig. [0(b), the 
electric field and equilibrium magnetization are along y, 
the displacement of the Fermi surface produces a non- 
equilibrium spin-orbit field (IB along x and the spin pre¬ 
cession around <5B produces a spin density aligned along 
z. The latter results in an additional torque that has 
a disorder-independent origin— This simple picture can 
be extended to more complex spin-orbit coupling situ¬ 
ation and only requires inversion symmetry breaking in 
the system. 


(a) 




FIG. 1. (Color online) (a) Fermi surface of a non-magnetic 
Rashba 2DEG: under the application of an external electric 
field (thick grey arrow), a non-equilibrium field <5B is pro¬ 
duced (thick blue arrow) that distorts the spin direction out 
of the plane (red arrows). After averaging, the spin density 
vanishes, (b) Fermi surface of a magnetic Rashba 2DEG in 
the strong ferromagnetic regime: in this case, since the spin 
directions (pink arrows) are initially mostly aligned along the 
magnetization m (black arrow), the resulting non-equilibrium 
spin density (red arrows) does not vanish and is aligned along 
z. 


III. TWO-DIMENSIONAL RASHBA 
FERROMAGNET 


We first apply this formalism to a ferromagnetic 2DEG 
in the presence of Rashba spin-orbit coupling.— This sys¬ 
tem is the prototypical free-electron model for SOTs in 
ultra thin ferromagnets embedded between two asym¬ 
metric interfaces! 7 ’ 8 ’ 10 Although the actual band struc¬ 
ture of magnetic bilayers such as Pt/Co is complex, re¬ 
cent first principle calculations indicate that this simple 
Rashba model qualitatively captures most of the relevant 
physics at these interfaces— This section is therefore de¬ 
veloped mostly for pedagogical purposes in order to make 
the dependence on various parameters explicit. The un- 
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perturbed Hamiltonian in Eq. m can be rewritten as: 


„ ^ 2^,2 

# 2 DEG = w — r - aa ■ (z x k) + J ex m • <x. (6) 

2m* 

where k = ( k x ,k y ,0) = k(cos(pk,simpk,0), a is the 
Rashba parameter and the magnetization direction is 
m = (cos(/?sin0,sin<^sin(9,cos(9). By diagonalizing 
Eq. ©, the eigenvalues and eigenvectors of itinerant elec¬ 
trons are 

£ k± = ^4 ± A k , 


2m* 


|k, +) — 


e I7k cos y 
sin 


;|k,-} = 



(7) 

- <p fe )sin0, 

(8) 

—e* 7k sin 
cos ‘y 

)(9) 


where we have cos \k = J ex cos0/A k and tan 7 k = 
ak cos tpk+J^ sm <p sm 6 ve i oc jty operator is given by 

OLK, Sin — Jex COS tp Sin C7 J S' O «/ 

v = ^Vk + x <t. Exploiting Eqs. ©-© in the weak 
exchange limit (akp J ex T), the nonequilibrium 
spin density reads 


ssintia = i w^ zxeE ) 

(5S i , nter = — (z x eE) 

1 8tt olE f v ; 

(m • z)eE 


inter 1 '^ex 


(5S‘ 2 r = 


4-7T CvEf 

and in the strong exchange limit (J ex akp T) 
1 am* 


^gintra _ 


2tt fi 2 r 


x [(z x eE) x m] 


1 ryrn * 

6S? tei = X [( z x eE ) X in] 

ro inter ^ OLTfl . . 

{Sj =-2i-r?r™ x ( * x eE) 


( 10 ) 

(ID 

( 12 ) 

(13) 

(14) 

(15) 


In summary, the SOT field defined as h = 2J ex SS/'yhN m , 
takes on the following form in the two limits: 


J ex <C akp : h = 


J PX am* 


2n"/N m h 3 T \ a 2 kp 


1 - 


r 2 


(z x eE) 


+ 


J 2 


2TT"/N m haEF 


(m • z)eE, 


Jex ^ akp . h — 


J ex am* 

'K^N m h 3 T V Je> 
am 


njN m h 3 


(16) 

1- y ) m x [(z x eE) x m] 

^ex J 

m x (z x eE) (17) 


field odd in magnetization direction. It lies perpendicular 
to the (m, z x E) plane in the strong exchange limit, see 
Eq. (1171) . This term is independent of the exchange J ex , 
in sharp contrast with the ISGE-induced SOT field, while 
in the weak exchange limit, Eq. (usd, it is second order 
in exchange and proportional to m z e E. The resulting 
intrinsic torque is in-plane in the strong exchange limit 
and even in magnetization direction. As will be seen in 
the next section, the parameters dependence displayed in 
Eqs. m and (fT71) is not restricted to the simple case of 
the Rashba model. Third, notice that in the strong ex¬ 
change limit the ratio of the anti-damping-like to the field 
like torque is ~ T/ J ex . This dependence is the inverse of 
what was found in Refs. □ and H3 in the diffusive limit 
and ignoring the interband scattering, where the ratio 
between the two torques is governed by J ex /r. A cor¬ 
rective anti-damping-like torque proportional to T s f/J ex 
is also obtained when considering a finite spin-flip relax¬ 
ation time T s f = fc/r s f . n i 14 

In the case of the anomalous Hall effect, related and 
better explored transport phenomenon, the intrinsic con¬ 
tribution dominates over the extrinsic contributions in 
the strong scattering limit i 39 ! 41 As a consequence, one 
is tempted to anticipate that the intrinsic contribution 
to the SOT discussed presently becomes important when 
strong momentum scattering is present (such as in disor¬ 
dered Pt/Co interfaces for example) and dominates over 
the corrections found in Refs. [ 11 ] and 0 in this limit. 
Nevertheless, these different contributions have been de¬ 
rived in different limits — i.e. strongkkll versus weak 
scattering (this work) — and should to be treated on 
equal footing for a rigorous comparison (e.g., see Ref. 
l4ll b Such a comprehensive model is beyond the scope of 
the present work. 


IV. DILUTE MAGNETIC SEMICONDUCTORS 
A. Method 

We now extend the previous results beyond the simple 
ferromagnetic 2DEG model with Rashba spin-orbit cou¬ 
pling. We consider a bulk three-dimensional DMS, such 
as (Ga,Mn)As, with a homogeneous magnetization. In 
order to model the SOT field of (Ga,Mn)As, we adopt 
a Hamiltonian including a mean-field exchange coupling 
between the hole spin (J) and the localized (d-electron) 
magnetic moment fjS a m of ionized Mn 2+ acceptor o 42 i 43 
and a four-band strained Kohn-Luttinger Hamiltonian. 
The total Hamiltonian of the DMS reads 


Three important facts ought to be pointed out. First, the 
extrinsic contributions (either intra- or interband) both 
give rise to an in-plane SOT field [even in magnetization 
direction, lying in the (m,z x E) plane]. The resulting 
extrinsic torque is then out-of-plane and odd in mag¬ 
netization direction. Second, the intrinsic contribution 
[second term in Eqs. (ITU and (1171) ] only produces a SOT 


J^DMS — T -^strain T Jpd-^MnSa^l ‘ J (1$) 

where J p d = 55 meV • nrn 3 is the antiferromagnetic cou¬ 
pling constant between hole and local moment spins for 
(Ga,Mn)As and S a =5/2 is the localized Mn spin. The 
hole spin operator is a 4 x 4 matrix^ The concentration 
of the ordered local Mn 2+ moments Nmu = 4x/a 3 is the 
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product of x that defines the doping by Mn 2+ ions and 
inverse volume per Ga atom (a is the GaAs lattice con¬ 
stant). The Kohn-Luttinger Hamiltonian in Eq. (fTSl) is 
expressed as^i 




- 

2 m 


7 lfc 2 I - 4:'Y3[kxky{J X , Jy} + C.p.} 

- 2 72 [(J| - \j 2 )k 2 x + c.p.] 


(19) 


This Hamiltonian applies close to the T point to centro- 
symmetric crystals with a diamond structure and strong 
spin-orbit coupling in the valence bands. The Luttinger 
parameters for GaAs are ( 71 , 72 , 73 ) = (6.98,2.06,2.93), 
I is the 4x4 unity matrix, J^,, J y , and J z , are the angular 
momentum matrices for spin |. They follow the relation 
{Jx, J y} = (J^Jy + JyJx)/ 2 , and c.p. denotes cyclic per¬ 
mutation of the preceding term. The first term denotes 
the kinetic energy of the holes. The second and third 
terms are associated with the spin-orbit coupling of the 
diamond crystal. In zinc-blende crystals, such as GaAs, 
bulk inversion asymmetry gives rise to the so-called cubic 
Dresselhaus spin-orbit coupling^ We neglect this term 
in the present study since there is no experimental in¬ 
dication that it contributes significantly to the SOT in 
(Ga,Mn)As. 

Hamiltonian Hdms should be understood as an ef¬ 
fective model attempting to describe the current-driven 
SOT in (Ga,Mn)As rather than the complete description 
of the electronic structure in this material. In a cubic 
diamond crystal, 72 7 ^ 73 . When 72 = 73 , dispersions 
following from Eq. m become spherically symmetric 
and when the spin-orbit coupling is removed completely 
(71 = 2 . 0,72 = 73 = 0), Eq. (fTf)l) reduces to a parabolic 
model. The impact of these three degrees of approxima¬ 
tion (parabolic model, spherical approximation and dia¬ 
mond crystal) on the SOT will be addressed in Section 
IIVB4I 

At this level of approximation the effective Hamilto¬ 
nian, Eq. in, does not break bulk inversion symmetry 
even though the actual crystal of the host GaAs does. In¬ 
deed, although the full model of GaAs contains additional 
terms that are odd in k (see Tab. 6.2 in Ref. flit ), it is 
experimentally established that the SOT in (Ga,Mn)As 
is sensitive to the strain. We therefore assume, in line 
with experiments^ that the key inversion-breaking term 
is proportional to the strain. The strain Hamiltonian is 
given by 

-^strain = C 4 [J x k x {e yy C zz ) 4” C.p.] 

T C§](Jx^y J y^x)Cxy T C.p. \ (20) 


where eu and e,; ; - ( i 7 ^ j) are the diagonal and non¬ 
diagonal elements of the strain tensor, respectively. We 
assume e xx = Cyy and e xy = Cyx • 

In Eq. (l20l) . the first term (oc C 4 ) originates from the 
lattice mismatch between the crystal structure of the sub¬ 
strate and the one of (Ga,Mn)As, and produces a spin- 
orbit coupling with Dresselhaus symmetry^! (oc e zz ). 


The second term (oc C 5 ) is the shear strain and possesses 
the symmetry of Rashba spin-orbit coupling^ (oc e xy ). 
Among the different terms linear in k and resulting 
from the inversion symmetry breaking (see Tab. C.5 in 
Ref. l46h . /^strain is the only one that acts in the mani¬ 
fold of heavy and light-hole states. It is worth pointing 
out that we consider here a large-enough system that 
allows us to disregard any effects arising from bound¬ 
aries and confinement. In the following, we assume 
C4 = C5 = 0.5 eV-n m 16 i 18 and consider the lattice mis¬ 
match strain (with e zz 7 ^ e yy = e xx — 0 and e xy = 0 ). 
Physical presence of the shear strain ( e xy 7 ^ 0) in unpat¬ 
terned (Ga,Mn)As samples is below the detection limit 4^ 
yet it has been introduced in previous theory studies 
to effectively model the in-plane uniaxial anisotropy ob¬ 
served in experiments^ Calculations with nonzero e xy 
are explicitly pointed out in the following. 

The SOT field h = 2J p d5S/"/h is evaluated once the 
energies E^ a and eigenfunctions IV'ka) implied by the 
Hamiltonian in Eq. (1181) are numerically calculated and 
the current-driven spin density dS is determined using 
Eq. @. In general, the SOT field can be decomposed as 

h = h m m + h\\e\\ + h±e± (21) 

where vectors e||,e_L have unit length, ej_||m x u, = 
e± x m, the subscript “so” has been removed for sim¬ 
plicity, and the direction of u (whose length is also set 
equal to one) should be chosen depending on the system. 
For example, we find u||z x E for the Rashba 2DEG. 
On the other hand, u = x in (Ga,Mn)As with growth 
strain (oc e zz ) as described by Eq. (TTSl) and current flow¬ 
ing along the [100] crystallographic direction. Our results 
presented below always assume ej_ pointing in the posi¬ 
tive z direction. 

In the following, we disregard the component of the 
SOT field which is parallel to the magnetization ( h m ) 
since it does not exert any torque on it. The two re¬ 
maining components in Eq. (1211) turn out to produce, in 
(Ga,Mn)As, the anti-damping-like SOT in the case of h± 
which is due to intrinsic interband mixing (of impurity- 
independent origin) and a combination of anti-damping¬ 
like and field-like extrinsic SOT in the case of h\\ which 
depends through T on the disorder strength. The angu¬ 
lar dependence of the two components, h\\ j_, reflects the 
details of the band structure as discussed in Sec. IIVB~4l 


B. Numerical Results 

For all the calculations presented in this section, the 
electric field E = 0.02 V/nrn is assumed to be applied 
along the x-axis and we varied the hole concentration 
between 0.3 nm —3 and 1 nrn -3 . This corresponds, re¬ 
spectively, to a Fermi energy between 200 and 450 meV. 
Except for Sec. IIVB 41 magnetization always lies along 
the y -axis (tp = 90°). 
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1. Intrinsic Versus Extrinsic Spin-Orbit Torques 

We first investigate the impact of impurity scattering 
on the intraband and interband contributions to the SOT 
fields. Figure [2] displays the SOT field as a function of 
the energy broadening T for different values of hole con¬ 
centrations. Although T is of the order of hundreds of 
meV in realistic (Ga,Mn)As, we choose T < 10 meV so 
as to be able to compare these results with the analytical 
ones obtained in Sec. m for the ferromagnetic Rashba 
2DEG which are valid in the small T limit. 

The intraband contribution to the SOT field, h|" tra , 
is inversely proportional to T for all hole densities as it 
is seen in Fig. HJa). This agrees with Eq. J3]) and also 
Eqs. m and (fldl) model for the ferromagnetic Rashba 
2DEG. No h™ tra component exists. On the other hand, 
the interband part (a ^ b ) of Eq. ([2]) contributes both 
to h || and h± which is shown in Figs. [2|b) and (c). The 
former is a correction to the intraband SOT field and it 
scales h™ ter oc T in the weak scattering limit. It tends to 
counteract the intraband contribution, as it is the case in 
the ferromagnetic Rashba 2DEG described by Eqs. m 
and dm The out-of-plane component h™ ter converges 
to a finite value when T vanishes, indicating the intrin¬ 
sic character of this part of the SOT field. These results 
are consistent with the analytical solutions obtained in 
Eqs. (fT3l) - (fl5l) in the ferromagnetic Rashba 2DEG and 
weak scattering limit. It is worth noticing that this de¬ 
pendence on spectral broadening holds over a wide range 
of T in the case of intraband contribution [see inset in 
Fig. Pfa)], while it breaks down already for T equal to 
few meV for the interband contributions. 


2. Ferromagnetic splitting 

The band structure of (Ga,Mn)As changes with the 
Mn doping that would, in the absence of the SOI, lead to 
a rigid mutual shift of the majority- and minority-spin 
bands. Such ferromagnetic splitting would be propor¬ 
tional to J ex = JpdNMnSa and we can distinguish two 
limiting situations in a system where the SOI is present: 
Eso T e x and Eso ^ T ex . In view of the analyti¬ 
cal results presented in Sec. cm it is meaningful to take 
Eso = ctk f in the Rashba 2D system. For each compo¬ 
nent of the non-equilibrium spin-density AS lntra , <5S 1 1 nter , 
(5S“ ter , there is a transition between different types of be¬ 
haviour in the two limits. For example, the out-of-plane 
component of the SOT field h changes from the oc Tg X be¬ 
haviour in the akp T ex limit implied by Eq. m into 
a J ex -independent behaviour in the opposite akp <§C T ex 
limit implied by Eq. (usd- We checked that this tran¬ 
sition occurs also in the numerical calculations across a 
range of J ex values. 

Contrary to the Rashba 2D system, the situation is 
more complicated in (Ga,Mn)As because of the addi¬ 
tional SOI terms in Eq. m- Due to their mutual 
competition, it is not obvious what should be taken for 



FIG. 2. (Color online) (a) Intraband and (b)-(c) interband 
contributions to the SOT field as a function of spectral broad¬ 
ening F for otherwise typical (Ga,Mn)As sample (doping con¬ 
centration x = 5%, lattice-mismatch strain e zz = —0.3%). 
Inset of panel (a) shows that /ij“ tra oc 1/T holds over a broad 
range of T. Only lattice-mismatch strain is considered, so 
that t xy = 0 in Eq. I|20l). 


the effective spin-orbit strength Eso- Looking at the 
J ex -clependence of the individual SOT field components 
in Fig. [3l we nevertheless recognize similarities to the 
Eso ^ T ex limit behaviour of the Rashba 2 D system. 
To some extent, this is a surprising finding since the dis¬ 
order broadening used for calculations in Fig. [3] is quite 
large (T = 50 meV), better corresponding to realistic 
(Ga,Mn)As samples but further away from the assump¬ 
tions used to derive the analytical results presented in 
Sec. IIII1 When J ex is small, both /i|| ltra and hj| lter are 
proportional to </ ex as seen in Eqs. (flOl) and (fill) , respec¬ 
tively. On the other hand, h™ teT oc Jg X in the bottom 
panel of Fig. [3] which is reminiscent of Eq. (fT31) . No 
similarities to the Rashba 2D system behaviour of the 
opposite limit (Eso <C T ex ) are found in our calculations 
for (Ga,Mn)As. 


3. Hole concentration 

We display in Fig. [I] the SOT field as a function of 
the hole density for different magnitudes of the lattice- 
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J ex< eV ) 


FIG. 3. (Color online) (a) Intraband and (b)-(c) interband 
SOT field as a function of exchange interaction J ex = JpdNMn- 
Varied values of J ex can be understood as a proxy to different 
Mn doping concentrations, e.g. x = 5% corresponds to J ex = 
0.06 eV, the spectral broadening is set to 50 meV and other 
parameters are the same as in Fig. [2] 


mismatch strain e zz . First of all, we notice that the SOT 
field components increase linearly with the strain. Sec¬ 
ond, increase of the hole concentration results in an in¬ 
crease in the in-plane SOT field h\\ approximatively fol¬ 
lowing a p 1 / 3 law, as shown in Figs. 0](a) and (b). This is 
consistent with Eq. (17) in Ref. [g in case of the intraband 
component. Interestingly, the in-plane interband SOT 
field h]j nter shows a similar tendency [Fig. 01(b)], while 

the out-of-plane interband SOT field h“ ter has a differ¬ 
ent dependence onp. This anti-damping-like SOT field in 
Fig • He) first increases with the hole concentration in the 
low hole density regime and later decreases towards a sat¬ 
urated value. This could be because of the competition of 
the different SOI types in (Ga,Mn)As as noticed by Kure- 
bayashi et al»22. Indeed, when the diamond-lattice spin- 
orbit coupling is absent (72 =73 = 0 ), the out-of-plane 
interband SOT field h^ ter increases with the hole con¬ 
centration following the same p 1 ^ 3 law as for the in-plane 
field [see inset of Fig. HJc)]. For a four-band Luttinger 
model that includes band warping (72 ^ 73 ), the compe¬ 
tition between the diamond spin-orbit coupling and the 


strain-induced spin-orbit coupling results in a reduction 
of h“ lter , as shown in Fig. 0]of Ref. (32j. The reason why 
the competition between the diamond spin-orbit coupling 
and the strain-induced spin-orbit coupling leads to the 
deviation from the analytical formula only in the case of 
the /i“ ter and not for / l ™ tra > lntel remains to be explored 
in detail. 

At this point, we remark that shear strain in Eq. JUD 
leads to /i™ tra comparable to values shown in Fig. [4ja) 
when the value of e xy is comparable to e zz used in 
Fig. 0J However, since the relevant values of e xy in unpat¬ 
terned epilayers are typically order-of-magnitude lower~ 
than those of e zz , we can typically expect an order-of- 
magnitude smaller /ij atra originating from the C' 5 -term in 
Eq. (P2U1) as compared to the C^-term. 



Hole density (nm 3 ) 


FIG. 4. (Color online) (a) Intraband and (b)-(c) interband 
SOT field as a function of hole concentration for different 
lattice-mismatch strain e zz . Inset in (c): interband SOT field 
in the parabolic model. The dashed lines in panel (a) are 
calculated using Eq. (17) in Ref. @ and follow ap 1 ^ 3 -law. Pa¬ 
rameters are the same as in Fig. 0]except for JpdNMn fixed to 
a value corresponding to Mn doping x = 5%. 


4- Impact of the Band Structure 

The total DMS Hamiltonian given in Eq. (fl8l) has both 
centro-symmetric and non-centro-symmetric components 
given by Eqs. © and mu, respectively. As discussed 

































in the previous section, the spin-orbit coupling of the 
centro-symmetric component of the Hamiltonian [i.e., the 
terms in Eq. m proportional to 72 and 73 ] affects also 
the SOT field, notably their dependence on the magne¬ 
tization direction [recall the definition of tp and 9 below 
Eq. ©]. Apart from the findings of Ref. [HI discussed 
above, it was shown in Ref. E that the shape of the 
Fermi surface has a strong impact on the angular depen¬ 
dence of the intraband SOT field /ij” tra . 

We now systematically explore the influence of the 
spin-orbit coupling of the diamond crystal on the dif¬ 
ferent components of the SOT field, i.e. foj” tra , /ij” ter 

and /i™ ter . The centro-symmetric component of the to¬ 
tal DMS Hamiltonian, Eq. m, accounts for the spin- 
orbit coupling through a set of the Luttinger parameters, 
71 , 2 , 3 - By tuning these three parameters, one can mod¬ 
ify the form of the centro-symmetric spin-orbit coupling. 
We model three distinct cases: (i) the parabolic approx¬ 
imation where no centro-symmetric spin-orbit coupling 
is present (71 = 2 . 0,72 = 73 = 0 ), (ii) the spherical ap¬ 
proximation where the centro-symmetric spin-orbit cou¬ 
pling is turned on but spherical symmetry is retained 
(72 = 73 = 2.5) and (iii) the diamond crystal where 
both cubic symmetry and centro-symmetric spin-orbit 
coupling are accounted for (72 ^ 73 ). This approach 
allows us to identify the role of the last two terms of 
Eq. (fTfll) on the SOT fields. In Fig. 0 we show the angu¬ 
lar dependence of the different contributions to the SOT 
field for the spin-orbit coupling induced by the lattice- 
mismatch strain in the context of models (i)-(iii). The 
magnetization lies in the (x,y) plane (9 = 7t/2) and its 
direction is given by the azimuthal angle ip. 

As expected from the symmetry of the C 4 term in 
Eq. (l20l) . the three components of the SOT field have 
a dependence of the form sin p in the parabolic model 
(□ symbols in Fig. 0). When diamond-lattice spin- 
orbit coupling is switched on but the spherical approx¬ 
imation is assumed, the interband SOT fields [A sym¬ 
bols in Figs. 0b) and (c)] deviate from this dependence, 
while the angular dependence of the intraband term re¬ 
mains unaffected [A symbols in Fig. 0a)]. Furthermore, 
the magnitudes of interband and intraband SOT fields 
strongly decrease. This is a manifestation of the com¬ 
petition between the strain-induced terms in Eq. (1201) 
with the centro-symmetric Luttinger spin-orbit terms in 
Eq. 

When the spherical approximation is lifted (72 ^ 73 ) 
electronic bands become warped, especially those of the 
heavy holes. This results into an increase of the interband 
SOT fields and an additional angular dependence shown 
by O symbols in Fig. 0 Microscopically, the latter ef¬ 
fect is caused by the distorted spin textures on the Fermi 
surface. The influence of the centro-symmetric spin-orbit 
field on the spin torque in GaMnAs has also been iden¬ 
tified by Haney et al<££ in DMS spin-valves. 



9 (deg) 

FIG. 5. (Color online) Intraband and interband SOT field as 
a function of the magnetization direction for different models 
labelled (i), (ii) and (iii) in the text. The red (O), blue (A) 
and black (□) data stand for the full four-band Luttinger 
model, its spherical approximation and the parabolic model, 
respectively. The parameters are the same as in Fig. [4] except 
for fixed p = 1.0 nm " 3 and t zz = —0.3%. 


V. SUMMARY AND CONCLUSION 

We have studied the intraband and interband SOT 
fields using Kubo formula, in the prototypical case of 
a ferromagnetic 2DEG with Rashba spin-orbit coupling, 
as well as in a three-dimensional DMS modelled by a 
kinetic-exchange Kohn-Luttinger Hamiltonian. For the 
latter, parameters pertaining to (Ga,Mn)As were used. 
In the limit of low doping concentration and weak ex¬ 
change coupling, we find similarities between the two 
systems, demonstrating that the general trends of the 
intrinsic and extrinsic SOT fields can be understood an¬ 
alytically using the Rashba 2DEG in the weak scatter¬ 
ing limit. Nevertheless, the numerical analysis of the 
three-dimensional DMS system also unravels the com¬ 
plex interplay between the different types of spin-orbit 
coupling (centro-symmetric and non-centro-symmetric) 
involved in realistic systems resulting in complex depen¬ 
dences of the SOT fields on the magnetization direction 
as well as significant differences from the Rashba 2DEG 
model. The contribution of interband mixing to the SOT 
presents an outstanding opportunity to explain the emer¬ 
gence of large anti-damping-like torques that cannot be 
readily attributed to spin Hall effect, offering an interest¬ 
ing platform to interpret recent puzzling results j 25 ’ 26 
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